The Magnetic-field Dependence of the Gaps in a Two-band Superconductor: A 
Point-contact Study of MgB2 Single Crystals 
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We present the results of directional point-contact measurements in MgB2 single crystals, in 
magnetic fields up to 9 T parallel to the c axis. By fitting the conductance curves of our point 
contacts - showing clear Andreev-reflection features - with a generalized Blonder-Tinkham-Klapwijk 
model, we were able to extract the values of the two gaps and A,r. The comparison of the 
resulting Ao-(i3) and A^(i3) curves to the theoretical predictions of the two-band model in dirty 
limit, recently appeared in literature, allows the first direct test of this model and gives a clear and 
quantitative proof that the tt band, even in the best single crystals, is in the moderate dirty limit. 

PACS numbers: 74.45.-l-c, 74.70.Ad, 74.25.0p 



More than two years after the discovery of supercon- 
ductivity in MgB2 Q the great experimental and the- 
oretical efforts of the entire scientific community have 
led to clarification of most features of this compound, 
mainly related to the existence of two band systems (cr 
and tt) and of the relevant gaps [3, 13 ■ Point-contact 
spectroscopy (PCS) has proved particularly useful in in- 
vestigating MgB2, since it allows measuring both the a- 
and TT-band gaps at the same time 0, |^ and determining 
the temperature dependency of these gaps with very high 
accuracy 0- However, some controversial points are still 
present. In particular, the effects of a magnetic field on 
the gap values and on the 7r-band contribution to the to- 
tal density of states (DOS) are still under debate, as well 
as the clean or dirty limit conditions of the two bands 
in the best samples and their influence on the measured 
physical properties. Both these topics can be addressed 
very well by performing accurate and directional PCS 
experiments in high-quality MgB2 single crystals in the 
presence of magnetic field. 

This paper presents some new results we obtained by 
PCS with magnetic fields up to 9 T apphed parallel to 
the c axis of MgB2 single crystals. These measurements 
gave us the magnetic-field dependence of the two gaps 
(Ajr and A^-) which results in good agreement with the 
theoretical predictions for the mixed state in a dirty two- 
band superconductor where the 7r-band diffusivity D-,: 
is a few times greater than the cr-band one, Do- 0. The 
consistency of this result is tested by comparing the zero- 
bias DOS calculated from our PCS data to the averaged 
zero-bias DOS predicted by the same model 0- 

The high-quality single crystals we used for our mea- 
surements were produced at ETH (Zurich) by using the 
growth technique described elsewhere . To control the 



direction of current injection, we chose only samples with 
regular shape, i.e. flat upper surface and sharp edges, 
even if these characteristics were fulfllled only by smaller 
samples (up to about 0.5x0.5x0.07 mm^). The point 
contacts were made by using either a small spot of silver 
conductive paint or a small piece of indium. The advan- 
tages of this technique with respect to the standard one 
(that employs a metallic tip pressed against the sample) 
have been widely explained elsewhere [3, ■ The normal- 
state resistance of our contacts ranges from 10 to 50 VI. 
Due to the rather large mean free path of our single crys- 
tals (€ 100 nm) |9|, they result to be in the ballistic 
regime |8|. 

In the experiments presented in this paper, the con- 
tacts were always placed on the narrow side of the crys- 
tals, so that, as already pointed out in previous papers 
0,01, the probability of electron injection was maximum 
along a direction parallel to the ab planes. Thus, in the 
following we will refer to our contacts as "a5-plane con- 
tacts" . The intensity of the magnetic field applied to the 
samples could range from up to 9 T and was affected 
by an uncertainty of about 1%, while the angle between 
the field and the ah planes was always = 90 ± 2° . 

The conductance curves {Al/dV vs. V) of our "spot" 
contacts were obtained by numerical differentiation of the 
I — V characteristics, and then normalized to the normal- 
state conductance. The upper curve in Fig^ (symbols) 
was obtained at T = 4.2 K in a Ag-paint contact. It 
has exactly the shape expected for a a6-plane contact 
13 , featuring symmetric conductance maxima at energies 
corresponding to A^- and A^ . This curve can be easily fit- 
ted by using a Blonder-Tinkham-Klapwijk (BTK) model 
generalized to the two-band case, where the total nor- 
malized conductance is given by cr = (1 — WT,)aa + Wt^Ut^ 
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FIG. 1: (a) Normalized conductance curves of a Ag-paint ab- 
plane contact, in increasing magnetic fields parallel to the c 
axis. Symbols: experimental data. Solid lines: best-fitting 
curves given by the two-band BTK model (see text). The 
curves are vertically shifted for clarity, (b) The same as in 
(a) but for a In-spot a6-plane contact on a different crystal. 



0, . Here, both and CTtt are calculated in the case 
of a spherical Fermi surface (FS), which is obviously not 
the case of MgB2. This problem cannot be overcome 
since, at present, there is no analytical theory for the 
Andreev reflection that takes into account the complex 
FS of MgBa. 

The fit has 7 free parameters: the gap amplitudes A^r 
and Act, the potential barrier coefficients and Za-, 
the lifetime broadening parameters T^r and To- and, fi- 
nally, the weight w^r- F is a phenomenological param- 
eter that usually accounts for pair-breaking effects (e.g. 
due to non-magnetic impurities) , broadening the conduc- 
tance curves and reducing their height with respect to 
the ideal ones. The best-fitting curve is shown in Fig^ 
as a solid line superimposed to the experimental data. 
The corresponding values of the fitting parameters are: 
A„ = 7.3 ± 0.1 meV, A^ = 2.8 ± 0.1 meV, = 0.944, 
Z^ = 0.484, = 3.3 meV, F^ = 1.46 meV, and finally 
Wtt — 0.75. The gaps are in very good agreement with 
the predictions of the two-band model (Eliashberg for- 
mulation) . The same happens for the weight , that 
is compatible with an injection cone of about 26° |5|. 



We can now try to fit in the same way the conductance 
curves measured in the presence of magnetic field. Due 
to the lack of an analytical theory of Andreev refiection 
for a two-band superconductor in magnetic field, this fit 
can be simply performed by using the lifetime broadening 
parameters Fo.,7r to mimic the effect of the magnetic field, 
as in Ref.Hol In this case, the total broadening parameter 
F is considered as the sum of an intrinsic Fj, due to the 
lifetime of quasiparticles, and an extrinsic Ff due to the 
effect of the magnetic field: F = Fj -f Ff. This approach 
assumes that the field breaks pairs in a way that can be 
represented by the quasiparticle lifetime, while its effects 
on the DOS are negligible in a first-order approximation. 

Unfortunately, at T = K the effects on the DOS 
curves of both the magnetic field and the actual shape of 
the FS are very relevant [Hill [11, so that our BTK- 
lifetime model fails in fitting them. We expect similar 
large effects to be also present in the Andreev-reflection 
conductance at T = K. Nevertheless, at low but fi- 
nite temperature (for example, T ~ 4 K) the thermally- 
broadened theoretical DOS curves become much more 
similar to those given by the standard BTK-lifctime 
model. To clarify this point, we tried to fit with this 
model some theoretical DOS curves at T = 4 K in the 
presence of magnetic field parallel to the c axis, i.e.: i) the 
local TT-band DOS curves calculated at various points 
of the vortex-lattice unit cell 6], ii) the two-band DOS 
calculated by taking into account the actual (simplified) 
shape of the FS of MgB2, made up of a distorted cylinder 
(for the a band) and a half torus (for the tt band) 
In both cases the lifetime-broadened BTK model proved 
very effective in fitting the theoretical curves in all the en- 
ergy range, apart from some deviations (< 30%) at very 
small energies {eV < 3 meV). The gap values determined 
from the fit are close to the theoretical ones used to gen- 
erate the curves, suggesting that Ff can mimic rather well 
the pair-breaking effect of the field at finite temperature. 

Thus, we can now proceed with the BTK fit of our 
experimental curves at S ^ 0. Up to about 1 T the 
fit is carried out by taking into account the contribution 
of the two bands, while above this field the fit can be 
obtained with no contribution of the tt band (i.e. by 
taking cr^ = 1). The best value of determined from 
the zero-field curve was kept constant in all the other fits. 
As in most junctions, also the values of Z^r and could 
be taken practically constant at the increase of the field, 
so that the actual free fitting parameters at i? / are 
Ao- , Att , Fo- and Ftt . The solid lines in Fig. QJi represent 
the curves that best fit the experimental data (circles). 
It is clear that the fit is very good, even if it cannot 
reproduce the asymmetry of some experimental curves. 

In Fig.^ the experimental conductances (circles) and 
the corresponding best-fit curves (solid lines) are shown 
for another junction obtained by pressing a very small In 
spot on the side surface of a MgB2 crystal. The best-fit 
parameters for the zero-field curve are: A^- = 7.0 ± 0.1 
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FIG. 2: (a) Symbols; magnetic field dependence of the gaps 
determined from the curves of FigCJi. The scale is on the 
lower and left axes. Lines: maximum gap amplitudes (cal- 
culated in Ref0 for B II c) as a function of the normalized 
magnetic field, when Da — 0.2_D.^. The scale is on the upper 
and right axes, (b) The same as in (a) but for the curves of 
FigE^- The theoretical curves are calculated for = Dt^. 



meV, = 3.1 ± 0.2 meV, = 0.51, = 0.41, = 
3.0 meV, F^r — 1.86 meV, and — 0.7. A comparison 
with the corresponding values of Fig^ shows that the 
big difference in shape between the curves in (a) and (b) 
is mainly due to the reduction of Zg- by a factor of two, 
while all the other parameters are practically unchanged. 
This is probably due to the different nature of the contact 
in the two cases (Ag paint in (a), indium in (b)). 

Figs-Et and report the field-dependence of the two 
gaps (solid circles) obtained by fitting the conductance 
curves of Fig. ^ and^jD, respectively. It is clearly seen 
that, for < 1 T, the 7r-band gap decreases at the in- 
crease of the field while the cr-band gap remains practi- 
cally unchanged. Above 1-1.2 T there is no longer trace 
of features related to in the conductance curves and, 
therefore, nothing can be concluded about this gap [i^ . 
The vanishing of the 7r-band contribution was observed 
in all the point contacts on MgB2 crystals we studied in 
magnetic field (more than 10) . This result is in complete 
agreement with many other experimental results present 
in literature 0,0, 0| and with the theoretical predic- 
tions of Ref. la, even if a very recent paper claims the 
possibility to determine the 7r-band gap up to about 5 T 
by PCS on MgBz films 0. 



In both Figs.EJ^ and[5j3, dashed lines represent the field 
dependence of the maximum pair potentials (A^g^^ and 
-^max) calculated at the boundary of the vortex-lattice 
unit cell for B || c 0. Even if our data rather represent 
an average of the conductance in different regions of the 
vortex lattice, the similarity between experimental data 
and theoretical curves is remarkable, provided that the 
overall shape of the curves is considered rather than the 
absolute gap values. In Fig|2t, the initial decrease of A^ 
is strikingly similar to that of AJ^^^^, in the case where 
Dg = 0.2Dt^. The behavior of A.^- in Figl^h is more lin- 
ear, and is here compared to the theoretical curve calcu- 
lated with Dg = Dt^. As far as Ao- is concerned, its field 
dependence agrees very well with the theoretical curve in 
both cases. In Figl^b, the theoretical curve well repro- 
duces the tendency of A^. to saturate at low fields. 

The internal consistency of the results reported so far 
can be checked by calculating the zero-bias DOS (ZBD) 
with the parameters we got from the fits, and compar- 
ing its field dependence with that reported in Ref^. To 
do this, i) we took the fitting parameters of each curve of 
Fig-H ii) we put ^^_cr = 20 and assumed a negligible con- 
tribution of the intrinsic lifetime broadening, i.e. we put 
T^''^ = 0; iii) we calculated with the lifetime-broadened 
BTK model the tunneling conductance curve at T = K, 
by using the values of A^-, A^, and F^''^ from the 
fit. The separation between intrinsic and field- induced 
pair-breaking effects is possible since the former are rea- 
sonably field-independent and can thus be determined 
from the fit of the zero-field curves. This procedure is 
somehow similar to that used in Ref^^ to evaluate the 
field-dependence of 7, which is in fact a thermodynamic 
probe of the low-energy DOS. 

In Figs.|3t and Oh the values of the ZBD determined 
from the data of Fig. ^ andQjD (solid circles) are com- 
pared to the total theoretical averaged DOS at eV=Q 
(solid and dashed curves) derived from the partial tt- and 
cr-band ZBD calculated in RefB In the case of Fig. |2t 
the agreement between the data derived from the ex- 
periments and the theoretical ZBD for — 0.2D-^ is 
impressive. On the other hand, the data of Fig. he 
between the total ZBD for D-^/Dcr = 5 and Dtj/D^ — 1, 
thus suggesting that an intermediate value of this ratio 
would allow a proper fit of the experimental data. Both 
these last results are fully consistent with those shown in 
Fig. 13 However, this agreement is not trivial since the 
quantities reported in Fig.|21and Fig. are almost inde- 
pendent. In fact, it can be easily shown that, for a given 
DOS curve, the value of the ZBD is little dependent on 
the gap amplitude, but strongly depends on the values 
of Ff' and F^. As a consequence, the consistency of the 
two different comparisons gives a strong evidence of the 
validity of our results. 

To summarize, these results indicate that: i) the life- 
time broadening can be used, as a first approximation, 
to mimic the effect of the magnetic field in the BTK fit 
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FIG. 3: (a) Symbols: magnetic field dependence of the zero- 
bias DOS (ZBD) calculated with Ao-.ir and F^''^ obtained by 
the fit of the curves of Fig[TJi, and taking Za = ^ 20. Solid 
line: field dependence of the theoretical ZBD for = 0.75, 
when Da ~ 0.20.^, 6j. (b) Symbols: ZBD calculated as in 
(a) from the curves of Fig^. Solid line: same as in (a) but 
with Wtt ~ 0.70. Dashed line: field dependence of the ZBD 
when Wrr — 0.70 and = D^r 0- 



gives a weighed average of A-^ over the FS. As far as the 
critical field i?c2 is concerned, we recently determined it 
by PCS for both B || c and B _L c HH. We always got 
values greater than those of Bc2 measured by bulk tech- 
niques (torque magnetometry, specific heat) and smaller 
than the surface critical field Bc3 HJl- The results will 
be widely discussed in a forthcoming paper, but we can 
anticipate that the high values of Bc2 measured by PCS 
are probably related to some non-trivial surface effects. 

In conclusion, in this paper we have presented the first 
detailed determination of the magnetic-field dependency 
of the gaps in a two-band superconductor by using the 
point-contact technique in high-quality MgB2 single crys- 
tals. The results allow the first direct test of the recent 
theory for the mixed state in a dirty two-band super- 
conductor. Moreover, they give clear and self-consistent 
evidence of the moderate dirty-limit conditions for the 
TT band and of the variability of these conditions from 
sample to sample. 

We are indebted to T. Dahm, O.V. Dolgov, A. A. 
Golubov and I.I. Mazin for discussions and suggestions. 
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02-02-17133), the Ministry of Science and Technologies 
of the Russian Federation (contract N. 40.012.1.1.1357) 
and the INTAS project N. 01-0617. This work was done 
within the project PRA "UMBRA" of INFM and the ASI 
contract N. I/R/109/02. 



of Andreev-reflection curves measured in MgB2 at finite 
temperature; ii) the field dependencies of A^-.tt and of the 
ZBD deduced from the experiments provide a clear and 
quantitative evidence that, even in the best single crys- 
tals, the TT band is in moderate dirty limit (Dj — 5-Do-), in 
complete agreement with recent STM [l5| and de Haas- 
van Alphen 1^ results on the same crystals, as well as 
with the model for electric transport in MgB2 jlBj and the 
theory of the mixed state in a dirty two-band supercon- 
ductor iii) different crystals, even if apparently pro- 
duced in the same way, may have different values of the 
diffusivity ratio D-^/D^, i.e. different degrees of cleanli- 
ness of the TT band. 

Some final comments are necessary on the zero-field 

lie . 

values of and on the upper critical field ini- 
plicitly reported in Fig|5| First, our PCS values of A^r 
(~ 2.7-2.8 meV)ll are systematically greater than those 
given by STM (~ 2.1 - 2.2 meV) Since the accuracy 
of both the techniques is out of discussion, this difference 
has probably a physical reason that is still under inves- 
tigation. Let us just remind here that the values of A^ 
are distributed between 1 and 3 meV over the FS |2Cll| . 
with a maximum at about 2 meV. Moreover, strong- 
coupling calculations |3| give an average value over the 
FS A^ = 2.7 meV. This might suggest that, while STM is 
more sensitive to the maximum of the distribution, PCS 
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